arXiv:1504.07804v2 [math.NT] 10Jun2015 


SUMS OF DIVISOR FUNCTIONS IN Fji] AND MATRIX 

INTEGRALS 


J.P. KEATING, B.. RODGERS, E. RODITTY-GERSHON AND Z. RUDNIGK 

Abstract. We study the mean square of sums of the fcth divisor func¬ 
tion dk (n) over short intervals and arithmetic progressions for the ratio¬ 
nal function field over a finite field of q elements. In the limit as g —^ oo 
we establish a relationship with a matrix integral over the unitary group. 
Evaluating this integral enables us to compute the mean square of the 
sums of dfc(n) in terms of a lattice point count. This lattice point count 
can in turn be calculated in terms of a certain piecewise polynomial 
function, which we analyse. Our results suggest general conjectures for 
the corresponding classical problems over the integers, which agree with 
the few cases where the answer is known. 


1. Introduction 

The goal of this paper is to study the mean square of sums of divisor 
functions over short intervals, for the rational function field over a finite 
field, and to use the results obtained to gain insight into the corresponding 
classical problem over the integers. 

1 . 1 . Classical theory. The fe-th divisor function dk{n) gives the number 
of ways of writing a (positive) integer as a product of k positive integers: 

(1.1) 4 (n) := #{(oi, ... ,afc) : n = ai • ... • Ofc, ai,...,> 1 } , 

the classical divisor function being d{n) = d2{n). 

Dirichlet’s divisor problem addresses the size of the remainder term A2(x) 
in partial sums of the divisor function; 

(1.2) A2(t) := X ^2(n) - ^(^logx -h (27 - 1)^ 

n<x 

where 7 is the Euler-Mascheroni constant. For the higher divisor functions 
one defines a remainder term Afc(x) similarly as the difference between the 
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partial sums Yln<x^k{n) and a smooth term xPk-i{logx) where Pk-i{u) is 
a certain polynomial of degree k — 1] see, for example, [38] Chapter XII. 

The mean square of A 2 (x) was computed by Cramer |9| for k = 2, and 
by Tong |39| for k > 3 (assuming the Riemann Hypothesis (RH) if A: > 4), 
to be 

2X 

(1.3) ^ j Akixfdx ^ CkX^~P 

X 

for a certain constant Ck- Heath-Brown m showed that Ak{x)/x^ has 
a limiting value distribution (for A; > 4 one needs to assume RH); it is 
non-Gaussian. 

1.2. The divisor function in short intervals. Let 

(1.4) Ak{x-, H) = Ak{x + H) - Ak{x) 

be the remainder term for sums of dk over short intervals [x,x + H], Our 
main concern is to understand its mean square. 

For relatively long intervals, Lester |29j proves an asymptotic (assuming 
RH for k > 3) similar to the result (|1.3p : 

2X 

(1.5) J (^Akix,H)ydx ^2ckX^-p 

The interesting range for us is that of shorter intervals: H < X^~k. For 
k = 2, Jutila [20], Coppola and Salerno [8], and Ivic mM show that, for 
X^ < H < the mean square of A 2 (x, H) is asymptotically equal to 

2X 

If/ \ 2 

(1.6) — J (^A 2 {x, H)j dx ~ iLF 3 (log—^) 

X 

for a certain cubic polynomial F^. In that regime, Lester and Yesha [30] 
showed that A 2 {x,H), normalized to have nnit mean-square using (II.6p . 
has a Gaussian value distribution, at least for a narrow range of H below 
the conjecture being that this should hold for X^ < H < for 

any e > 0. 

For k > 3, Milinovich and Turnage-Butterbaugh m p. 182] give an 
upper bound, assuming RH, of 

2X 

(1.7) ^ J (Afc(x, F)) ^dx « H{log , X^<H< X^-^ 

X 
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In concurrent work, Lester [29] shows that for k > 3, assuming the Lin- 
delof Hypothesis, if h{x) = 

( 1 . 8 )^ 

(21 o 1 

1 . o 2 -t: 1 


2X 


1 


a: 


2 1 9 2^~t — 1 

Ak{x,h{x))] dx ~ X*^-(logX)^ 


provided 1 — < d < 1 — -I, where 

(1.9) 

Tik + j)Y 1 


a-k = 


- LXJ 


i=o 


r(/c)j! J p 


p ^ i=o ^ ^ 


For k = 3 and ^ < <5 < |, the result is unconditional. 


1.3. A Conjecture. We did not find any conjecture in the literature for 
the order of growth of the mean-square of Afc(x; H) for small H. Based on 
Theorems oi and [LSI we believe the following: 

Conjecture 1.1. If0<6<l — ^ is fixed, then for H = , 

2X 

(1.10) ^ J (^Ak{x,Hfj dx ~ afcPfc((5)iL(log , X ^ oo 


where Ok is given by (HSl), and Vk{d) is a pieeewise polynomial function of 
6, of degree A:^ — 1, given by 

(1.11) P»{(S) = (1 - , 


Here 

( 1 . 12 ) 


7fc(c) = 


A:!G(1 + A:)2 


^ J 6ciwi +...+Wk)Y\i'>^i-'>^j)‘^d^w, 

[ 0 , 1 ]'= 


where dfix) = 6{x — c) is the delta distribution translated by c, and G is 
the Barnes G-function, so that for positive integers k, G{1 + k) = l\-2\- 
3\---{k-l)\. 

For 1 — < 5 < 1 — (jl.ll|l reduces to the simpler form 

rendering visible the compatibility of Conjecture 11.11 with Lester’s result 
m, which corresponds to taking H = with 6 in this range. Note that 

in m, the length of the interval h{x) = varies with x, and 

2- 1 

this slight difference in conventions is responsible for the factor of in 
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since the mean value of h{x) over [X, 2X] is 

2 ^ . 1 
1 f 2^~k - 1 . 

(1.14) - j h{x)dx = . 

X ^ 

As will be explained later, 7 fc(c) is a piecewise polynomial function of 
c, satisfying 7 a:(c) = 7 fc(A: — c), that relates to the asymptotics of a lat¬ 
tice counting problem iTheorem II.4p . This lattice counting problem itself 
emerges from the evaluation of a matrix integral over the unitary group. We 
also note that it is possible to write down conjectures for the lower order 
terms in the asymptotic expansion (II.ip : the right-hand side is, up to terms 
that are o(l), a polynomial in log A whose coefficients can be computed. 
This is explained in Section [5j 

1.4. Divisor functions in Fg[x]. We study the problem of the sum of 
divisor functions dk{f) over short intervals for Fq[x]. The divisor functions 
dk{f) for a monic polynomial / are defined in analogy to (II.ip and give 
the number of decompositions / = / 1/2 • • • /fc with /j monic. In particular 
d 2 = d is the classical divisor function. 

We denote by A4„ the set of monic polynomials of degree n. A “short 
interval” in Fg[x] is a set of the form 

(1.15) IiA;h) = {f:\\f-A\\<q^} 
where A G Ain has degree n, 0<h<n — 2 and the norm is 

( 1 . 16 ) ll/ll :=#IF,|*l/(/)=<(‘'**G 

The cardinality of such a short interval is 

(1.17) =: H . 

Set 

(1.18) MdM-^h):= dkif). 

fei{A-,h) 

The mean value is (c.f. U)) 

(ri!» ^ E = ^ E = ■ 

^ A&Mn ^ feMn ^ ' 

In analogy with (jl.2p and (jl.4p we set 

(1.20) Afc(A;h):=AA,,(A;/r)-g'^+i(”+;^^^) . 

We will show below lTheorem ll.2p that 

(1.21) = 


h> {1 — i)n - 1 

k 
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SO that Afc(j4; h) = 0 vanishes identically for h > {1 — ^)n — 1. The corre¬ 
sponding range over the integers is X^~k < H < X, where we have a bound 
of 0{X^~k'j for the mean square, see (ll.bp . 

Our principal result gives the mean square of Ak{A;h) (which is the 
variance of A/^j, [A] h)), in the limit q oo, in terms of a matrix integral. Let 
U he an N X N matrix. The secular coefficients Scj(JJ) are the coefficients 
of the characteristic polynomial of U : 

N 

(1.22) det(/ XxU) = Y^ Scj{U)x^ 

j=o 

Thus Sco{U) = 1, Sci(17) = trU, Scn{U) = detU. The secular coefficients 
are the elementary symmetric functions in the eigenvalues Ai,..., Aat of 17: 

(1.23) Sc,(17) = ^ Ail.A*, 

l<Zl < -<ir-<N 

and give the character of the exterior power representation on 

(1.24) Scj{U) =tTA^{U) 

It is well known that A-^ are distinct irreducible representations of the 
unitary group U{N), and hence one gets the mean values 

(1.25) j Scj{U)dU = 0, j = l,...,N 

U(N) 

and 

(1.26) J Scj{U)S^^dU = 

UiN) 


where the integrals are with respect to the Haar probability measure 

Define the matrix integrals over the group U(N) of Ax A unitary matrices 


(1.27) 


Ik{m-,N) := J 


Sc,Sc,,{U) 


dU . 


Then the variance 
(1.28) 


U(N) n+-+jk=ni 


1 


Var(A,J:=- |Afc(A;h)p 

^ AeMn 


satisfies 


Theorem 1.2. IfO<h< min(n — 5, (1 — ^)n — 2), then as q ^ oo 
(1.29) Var(A/'dJ = H • Ik{n;n- h-2) + . 
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In the remaining eases, 

(1.30) Var(AArfJ = h=[{l-^)n\-l 

and 

(1.31) Var(A/'rfJ = 0, [{1 - ^)n\ < h < n . 

In the case (1 — -^^)n < h + 2 < (1 — \)n, the matrix integral takes a simple 
form, c.f. Theorem 11.31 below. 


1.5. Matrix integrals. For [k — 1)N < m < kN, we obtain a simple 
formula for the matrix integral: 


Theorem 1.3. For {k — 1)N < m < kN, 

, , T- , fkN — m + k'^ 

(1-32) Ik{m;N)=\^ ^ 


We are also able to give a closed form, albeit more complicated, formula 
for the matrix integral for any range of the parameters, in terms of a lattice 
point count: 


Theorem 1.4. is equal to the count of lattice points x = G 

satisfying each of the relations 
(i) < x^p < N for all I <i,j <k 

(a) x^i'^ + x^ 2 ~^'^ + • • • + x^l'^ = kN — m, and 
(in) X G Ak, 

where Aj. is the eollection ofkxk matriees whose entries satisfy the following 
system of inequalities. 


™(1) 

Xi 

< 

^(2) 

Xi 

< ■ ■ 

■ < 

(k) 

VI 


VI 



VI 

™(i) 

< 


< ■ ■ 

■ < 

Jk) 

X2 

VI 


VI 



VI 

VI 


VI 



VI 

™(1) 

•^k 

< 

.^(2) 

< ■ ■ 

■ < 

(fc) 

4 


We note in passing that the above count of lattice points also may be 
interpreted as a count of plane partitions (see m, Section 7.20 for an in¬ 
troduction to the latter). 

For the standard divisor function (k = 2), if h < n/2 — 2 and n > 5 we 
thus find that as (? —>• oo. 


(1.33) 


4 

^ A&Mn 


2h + 6)(n 
"6 


2h + 7) 
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This is consistent with (jl.6|) . which leads ns to expect a cubic polynomial in 
§-/!. 

For the range (1 — ■j^)n < /i + 2 < (1 — ■^)re, (ll.32p gives 


(1.34) 




AeMr, 


k{n — h — 2) — n + k‘^ — 1 
k‘^ — 1 


1 , 


= HQf, 2 _i ((1 - -)n -{h + l 


where Q^ 2 _i{u) is a polynomial of degree k‘^ — 1, given by 

I — 1 

(1.35) 


Qfc2_i(u) . 


r(A;2) 


^ k^-i 

u'^ ^ + ... 


r(A:2) 


i- 1_ 

As this range corresponds to X '=-1 < H < X fc over the integers, the 
result (jl.34l) . (11.351) is comparable with Lester’s result (II.Sp (c.f. the remark 
after (ll.l3p L 

We use these results to model the situation over the integers for the range 
H < leading to Conjecture 11.11 To do so we derive asymptotics of 

Ikim; N) for m k, N: 

Theorem 1.5. Let c := m/N. Then for c E [0, k], 

(1.36) him- N) = ^uic)N^"-^ + OkiN^"-^), 
with jkic) defined by (|1.12l) . 

The matrix integral satisfies a functional equation /fc(m;iV) = hi^N — 
m-,N) (see Lemma l4.ip . from which it follows that 

(1.37) 7 fc(c) = jk{k - c) . 

It follows from an alternative analysis of /fe(m; N) that we also have 

Theorem 1.6. 

2 


(1.38) 


7 fc(c) = 






9k/ic - tj 


Q<l<c 

where gk,e{c — i) are (complicated) polynomials in c — £. 
and from this that 


Corollary 1.7. For a fixed k, 7 fc(c) is a piecewise polynomial function of 
c. Specifically, it is a fixed polynomial for r < c < r + 1 (r integer), and 
each time the value of c passes through an integer it becomes a different 
polynomial. 

For example, 

72(c) = ^ / (wi - ic - wi)f dwi = 

0<wi<l 
0<c—wi<l 


£_ 
3! ’ 


(2-c)3 


0 < C < 1 


1 < c < 2 
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and similarly 


73(c) 


0 < c < 1 
i(3-c)®, 2<c<3 


while for 1 < c < 2 we get 

73 (c) = ^(-2c®+24c^-252c®+1512c®-4830c'^+8568c^ 
8! V 


-8484c2+4392c-927 


1.6. Arithmetic progressions. A similar theory can be developed for 
sums of divisor functions along arithmetic progressions, see § [3l 


2. The divisor functions in short intervals 

Our first goal is to provide proofs for Theorem 11.21 and the other results 
on sums of dk in short intervals. 

2.1. An expression for the variance. To begin the proof of Theorem fL^ 
we express the variance of the short interval snms A/^j, in terms of snms of 
divisor functions, twisted by primitive even Dirichlet characters. Recall 
that a Dirichlet character is even if x(c/) = xif) for c G F®, and is odd 
otherwise. The number of even characters modulo is <h(r”“^)/(g — 

1) = (see e.g. [Ml §3.3]). We denote by — 1) 

the number of primitive even characters modnlo 
For a Dirichlet character y modulo set 

( 2 . 1 ) M{n-,dkx)-= X] dk{f)xif) ■ 

feMn 

Lemma 2.1. As q ^ 00 

(2.2) Var(V,.) = ^ ^ |A<(„; 4^)|2 + o( A) 

X mod T'^ ^ ^ 

X even primitive 

Proof. To compute the variance, we use |25l Lemma 5.4] which gives an 
expression for the variance of sums over short intervals of certain arithmetic 
functions a which are “even” (a{cf) = a{f) for c G ), multiplicative, and 
symmetric under the map f*{t) := in the sense that 

(2.3) a(r) = «(/), if/(0)/0. 

Since the divisor functions dk clearly satisfy all these conditions, we may use 
[25l Lemma 5.3] (compare [Ml §4.5]) to obtain 

TT ^ _ 

(2.4) var(A'rfj = - ^ 4(r™)4(r™) 

^ mi,1712=0 

^ 7^n-h\ dkx)M{m 2 ; dkx) 

X mod T'^-^ 

Xt^XO even 


X 
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To compute Ai{n;dkx)^ we introduce the generating function 


(2.5) M{m;dkx) 

m=0 


= 


Y. X{f)dk{f)u^^^^ = L{u,x)'^ 

f monic 

Hence Ai (n; dkx) is the coefficient of n” in L{u, x)^■ Now for even x 7 ^ XO) we 
write L{u,x) = (1 ~ u)P{u,x) and by the Riemann Hypothesis for curves, 
P{u,x) = ~ cijn) with the inverse zeros satisfying \aj\ < y/q. 

Hence we have an a-priori bound 

(2.6) |>f(m; 4 x)| <n,fc • 

Therefore in the sum (|2.4p . the terms with < 2n (i.e. (mi, m 2 ) 

(n,n)) will contribute 0{^q^~^) = 0{H/-^) (the coefficients 4 (T""'~™') = 
do not depend on q). Thus 


(2.7) Var(AA 4 ) = 


H 


4 


Y l-Ad(n;4x)P+ 0(;^) 


X mod T"- ^ 
X^XO even 


For the same reason, the non-primitive even characters, whose number is 
<C ^ev{T^~^)/q (see [Ml §3.3]), contribute 0{H/q) to the variance. Thus 
we are left with 

^ ^ = — _ i \ KA(-n-rl, ^ n( ^ ^ 


(2.8) Var(AAdJ = 


nify*lTn-h\ Y I-^(^;4X)P + 0( —) 

X even primitive 


□ 


2 . 2 . The sums M.{n]dkX)’ We need some information on M.{n]dkx) 

X even and primitive. By the Riemann Hypothesis (Weil’s theorem), for y 
even and primitive modulo we write 

(2.9) L{u,x) = {'i--u)det{I-uq^/'^e^) 
with G [/(n — h — 2) a unitary matrix of size n — h — 2. 

Lemma 2.2. For x even, primitive modulo T^~^: 

• If n < k{n — h — 2), that is h < {1 — ^)n — 2, then 

(2.10) M(n,x4) = (-l)V/' nS4(0x) + O(4-") 

jl+---+3k=n 

jl,---,jk<n-h-2 

• For k{n — h — 2) < n < k{n — h — 1), i.e. h = [(1 — ^)nj — 1 we get 

^ 

M{n,dkx) = 0{q~) 

• For n > k{n — h — l), that is h > (1 —^)n —1, we get M.{n,dkx) = 0- 


( 2 . 11 ) 
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Proof. For a primitive even character, the L-function is 
L{u, x) = (1 - u) det(/ - u^/qGy^) 

(ry 1 ^ ^ 2 

i=o 

To simplify notation in the calculations below, we write 
(2.13) N = n-h-2 , 


(2.14) Qj = (-l)-^g-^/^ Scj(0^), 0 < j < iV, a_i = 0 = qn+i , 

and set 

(2.15) bj := Uj — aj-i, j = 0,..., N + 1 
so that for x even, primitive 


(2.16) 

and 

(2.17) 


AT+l 

Hu, x) = X] 

j=0 


N+1 

Hu,x)^= E 

jl,..;jk=0 


Therefore the coefficient of tt"' in the expansion of L{u,x)^ for x even and 
primitive is 


(2.18) M{n, dkx) = E ^ji • • • • • 

jl+...+jk=n 

0<ji,...,jk<N+l 

Note that 

(2.19) bj = ttj + 0{qH ~)j 0 < j < N 
while 

(2.20) |6jv+i|<g^/^ 

Hence for an A:-tuple {ji,... ,jk) where one of the ji = N + 1 we have an 
upper bound 

( 2 . 21 ) . 


Thus if n > kN, and ji + ■ ■ ■ + jk = n, there is at least one index i so 
that ji = N + 1 and in that case 

_ 1 _ 

(2.22) \Min,dkx)\<-q~, kN < n < k{N + 1) . 
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For n < kN, there will always be an fc-tuple of 0 < ji,..., jfc < iV with 
ji + ■ ■ ■ + jk = n, and so for n < kN 
(2.23) 

M{n,dkx)= X] - 

jl+...+jk=n 

= (-!)>"/' Sc,,(0^) • ... • Sc,,(0^) + o[q^) . 

jl+---+jk=n 

0<jl,--;jk<N 

This concludes the proof. □ 


2.3. Proof of Theorem II.21 Inserting Lemma \n\ into (|2.2I) we find that 
for h<{l- \)n - 2, 


(2.24) 

Var(A/dJ 


H 


E 

X mod ^ 
even primitive 


jl + ---+jk=n 

0<ji,---,jk<n-h-2 


+ 0 



We now apply Katz’s equidistribution theorem for primitive even char¬ 
acters modulo |22j . which says that the corresponding Probenii 0^ are 
equidistributed in the projective unitary group PU{N — 2) if iV > 5, to 
replace the average over primitive even characters by a matrix integral over 
PU{n — h — 2), with an error of 0(l/y^). This gives 

(2.25) 


Var(A/'rfJ = H-Ik{n]n-h 



0 < /i < min(n — 5, (1 




which proves the main statement of our Theorem. 
In the remaining cases, Lemma 12.21 gives 


(2.26) Var(AA,J = 0(^), h=[{l-^)n\-l 

and 

(2.27) Var(A/'4) = 0, [(1 “ ■^)^J < ^ ^ • 

This concludes the proof of Theorem 11.21 □ 


3. The divisor function in arithmetic progressions 

3.1. Arithmetic progressions. We now turn to sums of divisor functions 
over arithmetic progressions. Set 

(3-1) 54 (A) = 54;X;q(^) = ^ 4(n) 

n<X 

n=A mod Q 
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For the standard divisor function {k = 2), it is known that individually, if 
Q < then 

(3.2) + 0(X>''3+»(i)) 

for some linear polynomial pq. This is apparently due to Selberg (unpub¬ 
lished). For recent work on asymptotics of sums of over arithmetic pro¬ 
gressions, see m and the literature cited therein. 

The variance Var(5d2;X;Q) of Sd2 has been studied by Motohashi [32] . 
Blomer [3] , Lau and Zhao [28] , the result being [28] (we assume Q prime for 
simplicity): 

i) If 1 < Q < XI/2 +e then 

(3.3) Var(5,2;A;Q) « + (^)'/'+^ 

ii) For <Q <X, 

(3.4) Var(5,2;X;Q) = ^P3(log ^) + 0((^)'5/6(logX)3) 

where ps is a polynomial of degree 3 with positive leading coefficient. See 
also the recent papers by Fouvry, Ganguli, Kowalski, Michel [H] and by 
Lester and Yesha m discussing higher moments. 

For k >3, Kowalski and Ricotta m considered smooth analogues of the 
divisor sums Sdf.-x-,Q{X), and among other things computed the varianc^il 
for < Y < 

We turn to Fg[x]. For Q E Fq[x] squarefree of degree at least 2, and A 
co-prime to Q, set 

(3.5) Sd,,n,Q{A) := 4(/) . 

f=A mod Q 


Our main result here concerns the variance 


(3.6) 


XaiQ{Sd^,n,Q) 


1 

HQ) 


Sdf,^n,Q{A) {<Sdi,,n,Q) 


A mod Q 
gcd(A,Q)=l 


2 


in the range n < /^(degQ ~ 1). 


Theorem 3.1. If Q is squarefree, and n < k{degQ — 1), then the variance 
is given by 


(3.7) 


. XaTQ{Sd,,n,Q) 
qn/\Q\ 


4 (n; deg (5 - 1) 


^The statement of |271 Theorem A], which deals with all moments, includes a term 
which is not small for the second moment; however the actual proof, see [271 equation 9.8 
and below] does give a good remainder. 
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In particular for the classical divisor function d = d 2 , we get a result 
consistent with (l.l.djl : 

Corollary 3.2. If Q is squarefree, degQ > 2 and n < 2(deg(5 — 1), then 
(3.8) 

Poh{n), n<degQ —1 

Pol 2 ,{ 2 {degQ - 1) - n), degQ <n< 2(degQ - 1) 

where Poh{x) = = (x + l)(x + 2)(x + 3)/6. 

As in the short interval case, we are led to a conjecture on the asymptotics 
of the variance over the integers. For simplicity, we stick with the case that 
the the modulus Q is prime: 

Conjecture 3.3. For Q prime, < X < as A —)• oo, 

Var(5dj^;X;Q) ~ 

where is given by (|1.9I) and 7 fc(c) is given by (I1.12p . 


lim 

q^oo 


VarQ(5d,,n,Q) 

q^/\Q\ 


3.2. Proof of Theorem 13.11 We start with the following expansion, using 
the orthogonality relation for Dirichlet characters to pick out an arithmetic 
progression [25l §4.1]: 

(3.9) 5d,,„,Q(A) =-^ ^ 4(/) + .^^ ^ ^(7l)7t^(n;4y) 

XT^XO 

(f,Q)=i 


where A4(n; dfcx), given by ( 12 . 11 ) . is the coefficient of in the expansion of 
L{u,x)^- Since L{u,x) is a polynomial of degree < degQ — 1 for y 7 ^ xo^ 
we see that Sd,^^n,Q is independent of A for n > fe(degQ — 1 ): 


(3.10) 


'Sdk,n,Qi-^) — 


HQ) 


E 

feMn 

(/.Q)=i 


4 (/) 


I (n+k—V\ 
1 fc-1 ) 

HQ) 


Thus for any n, the mean value (averaging over A coprime to Q) is 


(3.11) 


{•Sdk,n,Q) 


I fn+k—l\ 

1 k-1 J 

HQ) 


The interesting range is n < A:(deg Q — 1), which we assume from now on. 
To compute the variance, we use (|3.9p and the orthogonality relations for 
Dirichlet characters as in |24p 125] to hnd 

(3-12) Var(Sdk,n,Q) = \M{n]dkx)\^ ■ 

X^XO 

We first dispose of the contribution of even characters, whose number is 
‘hev(Q) = ^iQ)/{q — 1): As in (12.611 . we have an a-priori bound for x Xo 

(3.13) |A4(n,4x)| 4^2 ^ 
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Therefore the even characters contribute at most 


(3.14) 


<r 




r$ev(Q)g" < 


1 g" 


qm ) 


which is negligible relative to the main term that we find which is of order 
q'^/^{Q). The same argument bounds the contribution of odd non-primitive 
characters if Q is non-prime. Thus 


1 _ ( 7 ^ \ 

(3.15) Var(5rf,,„,Q) = ^^ ^ \Min;dkx)\^+o(^--^—) . 

^ X odd and primitive 


To handle the odd primitive characters x, we use the Riemann Hypothesis 
(Weil’s theorem) to write 

(3.16) L{u, x) = det(/ - , 

with the unitarized Frobenius 0^ E f7(deg(5 ~ !)• Hence for 2 < n < 
k{degQ- 1), 


(3.17) M{n;dkX) = (-!)>"/' Sc,,(0^) • ... • Sc,,(0^) . 

jl+---+jk=n 

0<Jl,...,ifc<degQ-l 


Inserting (I3.17P into (|3.15l) and using (|3.17l) and Katz’s equidistribution 
theorem |21] (here we require Q squarefree) we get for degQ ^ 2 and 2 < 
n < A:(deg Q — 1) 


(3.18) 

Var(5d,,n,Q) 

q-/\Q\ 


/ 


E 


Sc,,{U)...Scj,{U) 


^(degQ-l) 0<if!.tK<d^“gQ-l 
= 4(n,degQ - 1) , 


dU 


proving Theorem 13.11 

Note that If n < deg Q, then we of course do not need these powerful 
equidistribution results, since there is at most one f with deg f = n and / = 
A mod Q, which allows one to obtain the claim in an elementary manner. 


4. Matrix integral 

Our goal in this section is to evaluate the matrix integral (ll.27|) . We start 
by looking at the following products: 

(4.1) 

N N 

det(I - xU)’^ det(/ - yU*)’^ = ( Scj{U){-xy)HY^ Sci{U*){-yyy 

j=l i=l 

We will be interested in the expected value over the unitary group of the 
above. Due to the invariance of Haar measure of U{N) under multiplication 
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by unit scalars, we are left with only the diagonal terms, i.e. 

(4.2) f det{I-xU)’^det{I-yU*)’^dU = ^ Ik{m,N){xyr 

U{N) 0<m<kN 

This integral therefore serves as a generating series for the function Ik{m] N). 

Note that we may switch the sign of both x and y and retain the same right 
hand side. 

4.1. Evaluation in a certain range. We now give the proof of Theo¬ 
rem [T31 For the range m < N, we will apply the method of Diaconis- 
Gamburd m to obtain 

(4.3) 1 m<N 

When (k — 1)N < m < kN we have a functional equation which allows us 
to compute the integral in this range. 

4.1.1. The functional equation. 

Lemma 4.1. For 0 < m < kN, the following functional equation holds, 

(4.4) h{m-N) = h{kN -m-N). 

Proof. We use the functional equation of the characteristic polynomial of a 
unitary matrix 

(4.5) det(/ -|- xU) = x^ det([/) det(/ -|- x~^U*) 
which implies that 

(4.6) Scj{U) = det(17) ScN-j{U*) = det{U)ScN-j{U) 

Therefore 

Scj,{U) • • • Scj,{U)Sci,{U)---Sci,{U) = ScN-i,{U) • • • ScN-i,{U)ScN-n{U) • • • ScN-j,{U). 
We change variables 

(4.7) m' = kN-m, j'= N - ji, l[ = N-k 

and so obtain 

(4.8) 


E Sc,i(C/). 

.Scj,{U) 

2 

E Sc,;{(/)... Sc,-([/) 

il+-3k=m 




0<ji,...,jk<N 



0<j[,...,jl<N 


which implies (14.4p . 


□ 
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4.1.2. Review of Diaconis and Gamburd |12] . Let A = (aij) be an m x n 
matrix with non-negative integer entries. Let r* = mj be the sum of the 
entries in the Lth row, and cj = Y2i ^i,j be the sum of the entries in the j-th 
column. Set 

(4.9) row(A) = (ri,... ,rm), col(A) = (ci,...,c^) 

Let A = (Ai,..., Ar) G N'’ with Ai > A 2 > • • • > A^ be a partition of n, so 
n = Yhi Denote by m* = mi{\) the number of part of A equal to i, so an 
alternative notation is 

(4.10) A = (l”*i2™2...) 

Given two partitions /i = (^ 1 ,..., /im) and fi = (/ii,..., fin), denote by 
the number of m x n matrices A with non-negative integer entries so that 
row(A) = /X and col(A) = ft. For instance if ^ = (2,1,1) = (1^2^) and 
fi = (3,1) = (1^3^) then = 3 with the corresponding matrices A being 

2 0\ /2 0\ /I 1 

10 , 01, 10 
0 1/ \1 0/ \1 0 

We quote a result of Diaconis and Gamburd: 

Theorem 4.2. [12] Let ai,bj be non-negative integers, 1 < i, j < i. Then 
for max(X;^^i jaj, Ej=i jbj) < N, 

r ^ 

(4.11) j HiSc,iU)r Sc,{uf^dU = 

U{N) 

where /x = (l“i 2“2 . fx = (l'>i 2^2 _ 

4.1.3. Back to the variance calculation. There is a slight reformulation of 

Theorem 14.21 that will be useful to have stated. Let /x = {ji, and /x = 

(ji) ■■njk) be arrays of non-negative integers (we now impose no condition 
that they be weakly decreasing), and we generalize in the obvious 

manner, so that it is the count olkxk matrices A with non-negative integer 
entries such that row(A) = fx and co^A) = /x. Note that, by permuting rows 
and then columns of the matrix A, if the arrays n and v are rearrangements 
of each other, and likewise for /x and u, 

Thus Theorem l4.2l mav be reformulated as the statement that for max(^ ji, ^ ji) 

N, 

(4.12) J n Sc,^{U)'^JU) dU = 

U{N) * 

The reformulation is useful for us because in the proof that follows we will 
be working with arrays that are not ordered. 
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Proof of Theorem Al.Sl For m < N, note that in the definition (11.271) . the 
restriction that ji < N plays no role. Hence, 

2 


Ikim-N) : = 


UiN) 


Sc,,(17)... Sc,,(17) 

JlH- 'rjk=m 


dU. 


We may expand the square, and, because in the range of summation over j, 
we have Ji + • • • + jfc = m < A^, we may apply (|4.12p to see that the above 
expression is just 


E 




ii-i— 

h+---+3k=m 

But this sum is just the count of all /c x A; matrices comprised of non-negative 
integer entries with the total sum of the entries being m. This in turn is 
just the number of ways of writing ai 0^2 = m. Therefore, for this 

range of m < A^, /fc(m; N) is the binomial coefficient 

One way to see so is to note that it is the coefficient of x™ in 

1 


E 

ai>0 


X 


.CL\-\ -I"^fc2 _ 


(1 -x) 




Finally, to deal with the case {k — \)N < m < kN, we use the functional 
equation. Lemma l4.II □ 

4.2. Evaluation in other ranges. It was shown in the previous section 
how to evaluate /fe(m; N) in the ranges m < N and {k — 1)N < m < kN. 
Our goal here is to illustrate a general method for computing it in all other 
ranges. 

By (|4.2I) . we are looking to find the coefficient of x"* in the expansion of 


(4.13) 


Pf^(x) = J det(/ - U*xY det(/ - UydU. 

U(N) 


This can be calculated using the following Theorem: 

Theorem 4.3. [6], [7] Let A and B be finite collections of complex numbers. 
Then 


f det(/ - U*e-^) det(/ - Ue 

(M'l /3gS 


-^)dU 


U{N) 

(4.14) 


E ■ 

SCA 
TCB 
|T| = |S| 


3 -^(Eags “+Esct 


+ T-,T + S~) 
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Where 


S = A-S, T = B-T, S- = {-a, a G S}, T" = {-(3, /3 G T} 


and 


Z{A,B)= Hzia + P) 


with z{x) = . 

For example, we find 

P2{X) = 


a&A 


[1 + - (2 + iV) V+'V 


2 ,^1+Af 


(4.15) ' (l-x)4 

2(3 + 41V + N^)x^+‘^ - (2 + lV)^x^+3] 

Note that P 2 {x) satisfies x^'^P 2 (x)(l/x) = P 2 {x){x), which corresponds to 
the functional equation l 2 {m]N) = l2{2N — m;N). Evaluating the coeffi¬ 


cient of we recover 


(4.16) him; iV) = | Liv-l+si ^ ^ 


^2iV-m+3^ if AT < m < 21V. 


as proved in the previous section. 
Similarly 


P3ix) = 


1 


[1 - x3^+® + (3 + 1V)^(4 + 51V + lV^)(x' 


(1-x)^ 


2+7V 


— X 


Ar+7 


) + 


-(3 + lV)2(2 + Al)2(x2^+®-X^+^) + (3 + lV)2(10 + 71V + Af2)(^Af+4_^2Af+5)^ 

{3+Nf{N+4f{l/4){x^^+‘^-x^^^)+^{56+90N+51N^+12N^+N'^){x‘^^^^-x^^^)] 

Again P 3 (x) satisfies x^"^P 3 (x)(l/x) = P 3 (x)(x), corresponding to the func¬ 
tional equation 13 ( 1 x 1 ; N) = /3(31V — m;N). Hence 

’ (”^+®) if m < IV + 3 
Polyg(m) if IV -|- 3 < m < 21V — 3 


(4.17) 


h(m;N) = { 


(3^-™+®) if 21V - 3 < m < 3A^. 


Where Polyg(m) is a polynomial in m of degree 8, and is given by 


Polys ("i) = 
(3+A^)2(10+71V+1V2) 
(3 + Nf{4: + 5N + 


- (3 + Nf(N + 4f(l/4) 


m — N + 3 


+ 


m ^ +A_3^5g_^gQ^_^5^^2_^;^21v3-hlV^) 


m — IV -|- 5 




™-" + '5^_^(3 + Jvrt2 + iV)= 


m — IV -|- 7' 


This method obviously extends to larger values of k, but in practice is 
effective when k is relatively small. 
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4.3. Large N asymptotics: A symmetric function theory approach. 

In this subsection, we give a proof of Theorem 11.51 determining the asymp¬ 
totic behavior of N) when m and N grow in ratio to one another. We 
begin however with a proof of Theorem ll.41 the characterization of A") 
in terms of a count of lattice points. It is then in part by estimating this 
lattice count that we obtain the coefficient 7 fc(c) in Theorem 11.51 

4.3.1. Some preliminaries from symmetrie funetion theory. The proof below 
of Theorems 11.41 and 11.51 requires some knowledge from symmetric function 
theory. In order to make our presentation self-contained, in this section we 
recall for the reader a few concepts that will be necessary. In particular 
Schur functions, defined below, will play a key role. The reader already 
familiar with this material may skip ahead to the next subsection. (Standard 
references for this material include mmiEZ]; for readers with a background 
in analytic number theory, m is perhaps the quickest general introduction.) 

Recall (from ITTT2]) . a partition A is a sequence (Ai,...,Afc) of positive 
integers satisfying Ai > A 2 > • • • > A^. The length £{X) of such a partition 
is defined by £{X) := k. If 1 appears among the numbers Ai,..., A^ a total of 
mi times, 2 appears m 2 times, and so on, we also write A = ...). 

A Young diagram is a collection of boxes arranged in left-justified rows, 
with a weakly decreasing number of boxes in each row. The partition 
(Ai,..., Afc) corresponds to a Young diagram with Ai boxes in the first row, 
A 2 boxes in the second, and so on to A^ boxes in the fcth row. For instance, 
the partition (6,4, 3,1) corresponds to the Young diagram 



For A a partition, a semistandard Young tableau (SSYT) of shape A is 
an array T = (Tij)i<i<£(A),i<j<Ai of positive integers such that Tij < Tij+i 
and Tij < Tj+ij. It is common to write SSYTs in a Young diagram, as for 
example 


1 

1 

2 

3 

3 

7 

2 

3 

3 

4 



4 

4 

6 




7 







This is a SSYT of shape (6,4, 3,1). Note that the condition Tij < Tij+i 
translates to the array T weakly increasing in every row and Tij < Ti+i,j 
to strictly increasing in every column. 

T has type a = ( 01 , 02 ,...) if T has Oj = ai{T) parts equal to i. The 
SSYT above has type (2,2,4,3,0,1, 2). It is common to use the notational 
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abbreviation 


T _ “i('r) “2(T) 


•Xj - tX/ 1 iX/ 


SO for the example SSYT above, 

T 

X = 

We finally come to the combinatorial definition of Schur functions. 


T 2 2 4 3 
X = XiX2X^X4^XqX7. 


Definition 4.4. For a partition X, the Schur function in the variables 
xi,...,Xr indexed by X is a multivariable polynomial defined by 

sx{xi,...,Xr) ■■= 

T 

where the sum is over all SSYTs T whose entries belong to the set {1, ...,r} 
(i.e. afiT) = 0 for i > r). 

For example, the SSYTs of shape (2,1) whose entries belong to the set 
{1, 2,3} are 


1 

1 

1 

2 

1 

CO 

1 

1 

1 

2 

1 

CO 

2 

2 

2 

CO 

2 


2 


2 


3 


3 


3 


CO 


3 



and so 

S{2,1){xi,X 2,X3) = xfx2 + XiX^ + xfxs + + X^X^ + X 2 x| + 2 X 1 X 2 X 3 . 


4.3.2. A proof of Theorems \1.4\ and \l.^ 

Proof of Theorem \1.4\ Our starting point is again equation (14.21) . which in 
this case we evaluate using a result of Bump and Gamburd [5l Prop. 4]: 

Theorem 4.5. [5] Let ai, ...,aL+L' be complex numbers. Then, 

n det(l + af^U-^) n det(l + aL+eU) dU = _ 

e=i t=i '"“a 

Here is a Schur function indexed by the partition (N^). 

By specializing this Theorem, we see that, 

(4.18) 

det(l + all)^ det(l + dU = s ,..., /3,..., /3). 

U{N) k terms k terms 

Expanding the Schur function as a polynomial and labeling the coefficients, 
we have 

k terms k terms 

By comparison with (14.211 . we see that Ikinx] N) = Cmm- 


I 

U{N) 
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From the combinatorial definition of Schur functions lDefinition l4.4l abovel. 
we see that Cmm is the number of semistandard Young tableau (SSYT) T 
such that if, as before, ai denotes the number of i’s in T, 


flfc+i + • • • + a2k — rn, 

and Oj = 0 for i > 2k. 

We parametrize such tableaux T by letting = yi^\T) be the rightmost 
position of the entry s in row r; if s does not occur in row r, inductively 
define yr^^ = yl^ with yl^^ = 0 if the entry 1 does not occur in row r. So, 
for instance, in the SSYT T on the partition (7^) with entries ranging from 
1 to 4 given by 


T = 


1 

1 

1 

2 

CO 

CO 

CO 

2 

2 

2 

4 

4 

4 

4 


we have 



(3) 

Vi 

(3) 

2/2 


7 

3 


Note that here y^^ = 0 and y® = y[^'^ = y^'^ = 7. That these entries should 
take these values is necessarily the case; if the 2^^*^ row began with 1, then 
T could not be made to be strictly increasing in columns, and for the same 
reason the 1®* row may not end with 4. 

Moreover, note that because rows increase weakly, 

(4.19) 

and because columns increase strongly, 

(4.20) 

With these restrictions (I4.19P and f|4.20p in place, there is a bijection between 
arrays 


(y? 

0 


( 2 ) 

Vi 

( 2 ) 

2/2 


(fc) 

2/i 

(fc) 

2/2 


N 


N 


n\ 

N 


VO ^ vr vr’ • • • vr ’ 

with yi^^ € [0,and SSYT of {N^) with entries ranging from 1 to 2k. 
It is easy to see that those SSYT for which afc+i + - • ■+a 2 k = m correspond 

to those arrays in which {N — y[^'^) + {N — y^'^) H-+ (A^ — y^^^) = m. By 

re-indexing xi^'^ = we obtain the proposition. □ 


With Theorem 11.41 in hand, getting an expression for 7 fe(c) in Theorem 
fT5l as we will see, is a more or less standard argument in counting lattice 
points. On the other hand, in order to simplify the expression we get to 
(|1.12p . it will be useful to have done the following computation beforehand. 
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Lemma 4.6. As usual, define the Vandermonde determinant by 

A{wi,W 2, Wk) := - Wj), 

i>j 

and for fd G satisfying /3i < /32 < • • • < (3k+i, define 

/(/3) = {a € : /3i < ai < /32 < 02 < • • • < afc < Pk+i}- 

Then 

(4.21) j A{ai,a2,...,ak)d'^a = 

«e/(/3) 

Proof. Because of the well known identity A(i(;) = det(tt)J)“^), we see that 
the left hand side of (I4.2ip is just 




/ 1 1 

... 1 \ 




«! 02 

Ofc 


/ det 



d^a 

aeiih) 


a^- 

' ••• 



/ 

h - fii 

^3 - h 

Pk+l — fik \ 

= det 

(/3 

1 - /3?)/2 

iPi - /3|)/2 

•• (/?.%!-/?D/2 





•• iPk+i-Pk)/kJ 


by integrating one variable at a time and using multilinearity. But again 
applying multilinearity (twice), we see that this is just 

P 2 + £2 ■ ■ ■ Pk+ek\ 

Phe, ••• 

• . • ? 

••• 

where |e| is the number of i such that £i = 1. Clearly the determinant in the 
summand will be 0 unless s is one of the k + 1 possibilities: (1,1,1,..., 1), 
(0,1,1,...,!), (0,0,1,...,!), ..., (0,0,0,..., 0). Thus the sum above is just a 
Laplace expansion of 



/I 

1 • 

1 \ 

— det 

fii 

h • 

• fik+l 

kl 


/3| • 



1 


^ (_l)fc-Ndet 

£G{0,1}'' 






as claimed. 


□ 
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Proof of Theorem We demonstrate first that (jl.36jl of Theorem ll.5l holds 
with 7 fc(c) given by 

(4.22) lk{c) = j 5c{uf'' ^ - 'ru^^"*)lAk{u)d^‘^u, 

[ 0 , 1 ]'“^ 

where is the indicator function of the set (defined in the statement 
of Theorem [131) • 

The truth of this should come as no surprise; we have just approximated 
a lattice count with a continuous approximation. Later we show that this 
integral is equal to the right hand side of ()1.12p . 

Our proof of this first part is standard. For notational reasons let S = 
{{hj) : 1 < LJ < ^ : (bi) 7^ (1)^)}; let 14 be the convex region 
contained in ^ defined by the following 

system of inequalities: 

(i) 0 < uP < 1, for all {i,j) G S, 

(ii) For uP := c — + • • • + uP), we have 0 < uP < 1, and 

(hi) The matrix {uP)i<ij<k lies in the set Ak- 

This region is convex because it is the intersection of half planes. Note 
moreover that for all c € [0, fe], the region 14 is contained in [0,1]^ and 
therefore contained in a closed ball of radius y/k'^ — 1. 

Theorems O and Lemma l4. II show that 

(4.23) 4(m; N) = n {N • K)), 

where N ■Vc = {Nx : x G W} is the dilate of W by a factor of N. 

We will need to reference the well known principle that a count of lattice 
points in a region can be approximated by the volume of the region (at 
least in ordinary circumstances). A result of the sort we quote below dates 
back to Davenport mm-, the clean formulation we have cited here may 
be found in [Ml Section 2]. 

Theorem 4.7. If S C is a eonvex region contained in a closed ball of 
radius p, then 

(4.24) #(5 n Z^) = vokiS) + 
where the implicit constant depends only on 1. 

Applying ()4.24D . with £ = — 1, we see 

4(m; iV) = vol(fV • W) + ©^(iV^'-"). 

Yet clearly 

vol(iV • K) = J Sc{uP+ --- + uP)lA,{u)d’^\, 

[ 0 , 1 ]'“^ 


which implies ()1.36|1 . with 7 fc(c) given by ()4.22|1 . 
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It remains to show that this integral can be reduced to the expression 
defined in (jl.l2p . Here we make use of Lemma [4.61 We have, by applying it 
inductively, 


7fc(c) = j Sciu[^^ + 
[0,1]'=^ 




^ / 


4^’ < < 


VI 

< 


VI 


VI 

VI 


H-h • 1 


“4^ 

VI 





VI 


VI 



< 

(fe-1) 

“Li 

< 

{k) 

<-l 




VI 


VI 



< 


< 

(fc) 

“fc 


< 





VI 







< 





VI 









^ “fe 


VI 

k-l 

VI 

(fc-1) 




1 ! 


1 ! 


VI 

,(fc) 




[ 0 , 1 ]' 


^ H-h < • • • < 


l!-2 !---(A:-l)! 


l!-2!---(A:-l)! 


d^u 


dciwi + • • • + Wk)^{w)‘^ d^w, 


k\G{l + ky^ 


[ 0 , 1 ]* 


with the last step following from symmetry. 


□ 


We note for the reader familiar with Gelfand-Tsetlin patterns that what 
we have done in these last few steps is to compute the volume of what is 
called a Gelfand-Tsetlin polytope. A computation of this volume has ap¬ 
peared before in the literature (see [2] for a proof using representation theory, 
or |33] for a proof using the Harish-Chandra-Itzykson-Zuber integral), but 
the elementary proof we give here based on Lemma 14.61 seems to be new. 


4.3.3. Ehrhart theory. Theorem 11.41 also allows us to say something about 
the algebraic character of the quantities we have been discussing. 
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Corollary 4.8. Let c = p/q be fixed rational number and k be a fixed integer. 
If N is a multiple of q, then Ik{cN, N) = where Pc^k is a polynomial 

of degree k'^ — 1. 


Proof. This corollary follows from an application of a theorem of Ehrhart 

m- 

Theorem 4.9. If E is a convex lattice polytope in (that is, a polytope 
whose vertices are all integer coordinates), then there is a polynomial P of 
degree n, such that for all I G N>o, 

= P{i). 


Returning to the corollary at hand, we have from (I4.23p . when N = qi, 
Ik{cN;N)=#{Z'^"-^n{£-[q-Vc])). 

But then it is straightforward to verify that qVc = qVp/q is a convex lattice 

polytope in so that Ik{cN‘,N) is a polynomial in i and therefore in 

N. □ 

4.4. Large N asymptotic: the complex analysis approach. In this 
subsection we prove Theorem 11.61 The approach we take is based on the 
following expression proved in [6] (Lemma 2.1): 


Theorem 4.10. [6] Letai,l3j be complex numbers. Then, 


U{N) 


]Jdet(/ - n “ Ue-^i)dU 

i=i 


i=l 


qN {ai-\ -h«r) 

(27ri)2^(r!)2 

A(zi,...,Z2r)^ 


^ — - \-Z2r^ 


n (1 


l<i<r 

r+l<q<2r 


-dzi ■ ■ ■ dZ2r 


Where A{zi ,..., Z 2 r) = — zf) is the Vandermonde determinant, and 

the contour integrals enclose the variables 

From the definition ()4.13p . we have 
(4.25) 

Pk{x) = 

(_l)fc^fc7V 


{27ri)‘^^ {k\fi 

A{zi, . . . , Z2k f 

where a = logx. 


Af(2fc^iH- \-Z2k) 


n (1 


l<l<k 

k+l<q<2k 


dzi-■ ■ dZ2k 
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We set m = cN and will consider when 0 < c < A: is fixed and iV —>■ oo 
in such a way that cN is an integer. We will then need to compute the 
coefficient of in Pk{x). 

We first shrink the contour in (I4.10j) into small circles centered at 0 and 
a. This leads to a sum of 2^^ multiple integrals, each surrounding either 0 
or a; c.f. the calculation in [23]. Taking into account symmetries between 
the variables and counting the number of ways of picking i of the first k 
contours to surround a, and k — i oi the second k contours to surround a, 
we find 

(4-26) Pk{x) = Pk,i{x) 

£=0 ^ ^ 


where Pk/{x) is the integral with contours zi,... ,Zi, Zk+e+i, ■ ■ ■, Z 2 k along 
small circles surrounding a = logx and zi^i,..., Zk+e along small circles 
surrounding 0. The remaining integrals where there are different numbers 
of contours surrounding a and 0 do not contribute, as proved in the following 
Lemma 14.111 which we prove below. 

Next we change variables 


Vj 

where 

^ ^ f 1 , j = I,. oi j = k + i + 1,... ,2k 

~ \o, i+l<j <k + £ 

This gives that the integrand of Pk,i{x) is, up to terms of order 1/N smaller, 


g-K+i+...+^2fc) n a 2 n {^^?dvi...dv: 


'J2k 


(4.27) 


„-N{k-t) 


i<j 


i<j 


J\f2k 


_ 2k 


n N 

t<k<g t<k<g 


i=i 


The number of pairs i < j with 7 ^ ej is /c^, hence n i<j a? = and 

the number of pairs i < j with = ej is (^ 2 ^) — k'^ = k'^ — k, so that 


n 

Kj 


,Vi - V. 


j \2 


1 


N 


]Sl2(k^-k) 


n • 


i<j 

ei=ei 


The number of pairs (t, q) with l<t<k<q<2k and et 
hence 


n 

t<k<q 

^t—^q 


Vq - Vt 

N 


1 

]\[2l{k-l) 


n 

t<k<q 

et=eq 


Vt) . 


eq is 2i{k — i), 
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Therefore (j4.27p is equal to 

Kj i 

(_l)fc^-Af(fc-^)jY 2 <;(fc-€)_f 2 Z!£_. 

0(1- ) 0 (^5 - ^t) 

t<k<q t<k<q 

et/cg ^t=^q 

In the denominator, we rewrite the expression nt<A:<g(l — n ) by 

noting that is x if = 1 , et = 0, which happens when t = £ + 1,... , fe 

and q = k + I + 1,..., 2k, and it equals x~^ if = 0 and et = 1 , which 
happens when t = 1,... ,i and q = k + 1,... k + i. Thus 

k 2k i k-\-i 

t<k<q t=£-\-l q=k-\-£-\-l t=l q=k-\-l 

I k+e 

=i-iYx-^" n n 

l<£<fc t=l q=k-\-l 

k+l<q<2k 

Multiplying by the common pre-factor of ^ — gives that, up to a term 

of order 1/N smaller. 


(4.28) 

Pk,i{x) ~ (-1) 


^y.t{N+i)j^2e{k-i) ^ 

(A:!)^ (27ri)2*^ 


n (1 




l<t<k 

k+l<q<2k 

<^t¥=<^q 


I k+l 

n n 

t=l q=k-\-l 


'|g-(Dfe + l + ...+ 1 ) 2 fc) 


2k 




dvi 


l<t<£, k+l+l<q<2k 
or 

£-\-l<t<k, k-\-l<q<k-\-£ 


or 

k-\-i-\-l<i<j<2k 

or 

i-\-l<i<j<k 

or 

k-\-l<i<j<k+i 


i=i 


We need to pick out the coefficient of x’^^ in Pk/{x). (This coefficient 
is automatically 0 if £{N + t) > cN, so we need only consider £ < c.) We 
therefore need to find the coefficient of 


l<t<k 

k+l<q<2k 


(4.29) 

















28 


J.P. KEATING, B.. RODGERS, E. RODITTY-GERSHON AND Z. RUDNIGK 


We can expand the above to get 

OO 

(4.30) Y 

m =0 6i + ...+6^2_,_(fc_^)2=m 
bi>0 

If we consider the pre-factor of 04=1 nij=fc+i e ^ ^ then the required coef¬ 
ficient is 

(4.31) trSym("-^)^exp(^y) 

where V := diag(ug — vt) for q and t such that l<t<A:, k + 1 < q < 2k 
and et / Cg. Next, we use Lemma [4.121 proved below, to deduce that the 
expression (14.311) is 

(4.32) ((c - ^)iV)^'-2^(^-^)-i J^((c - £)v) 
with 

(4.33) Jeivi,...,V2k)= j 

^tq>0 

where the + {k — i'f' variables x^q have indices l<t<k,k + l<q<2k 
with et 7 ^ Cq, that is either l<f<£, /c-|-l<g<fc-|-£orl'-|-l<f<A:, 
k + i + l<q< 2k. 

Since Pk/{x) also has a factor of we get a total contribution of 

— i) where 

(4.34) g^c -l)= t)v) 

2k , 

^-{Vk + l+- + V2k) U i^y. _ y.-j2 ^ 

l<t<£, k+e+l<q<2k i<i<j<i j=l 

or or 

£-\-l<t<k, k-\-l<q<k-\-£ k+£-\-l<i<j<2k 

or 

£-\-l<i<j<k 

or 

k+l<i<j <k-\-£ 

The prefactor Qk/ic—i) depends polynomially on c—i, because to compute 
it we need to compute derivatives of Je{{c — i)v) at u = 0, which are clearly 
polynomial in (c — i). 

Summing these over 0 < £ < c gives an expression of the form 'Ykic)N^ 
where 

(4.35) 7t(c)= E 

0 <£<c ^ ^ 

as was to be proved. 

It remains now to prove the two lemmas we have used. This we do in the 
following subsections. 
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4.4.1. Vanishing of an integral. Denote by aei,..., ae 2 k) the integral 
Pk{x) over the circular contours centered in aei when €i can be either zero 
or one. 


Lemma 4.11. Let the number of ei which are equal to 1 and the number 
which are equal to 0 be different. Then the integral Pk{x;aei,... ,ae 2 k) is 
identically zero. 


Proof. We consider the case in which there are more zeros then ones. The 
case in which there are more ones then zeros, can be deduced in the same 
way. We can choose (without loss of generality) ei,..., Ck+i to be zero. 
Denote 


(4.36) 

G{zi ,..., Zk+i) := e-^(^'=+i+-+^ 2 D 

l<Z<fc 

k+l<q<2k 


-1 


A{zi, . . .,Z 2 k) 


ntlip-armk+2i^^) 


2k 
i=k+2\ 


This function is analytic around zero. The poles that arise when Zq = zi 
cancel with the Vandermonde determinant. Next, we use the residue theorem 
in order to compute the integral. Consider the Vandermonde determinant 
expansion: 

2k 

A{zi,...,Z 2 k)= Sgn(cj)(]J( 2 :i)'"h)-i) 

o-eS2k *=i 

By the residue theorem we need to show that the coefficient of 

in the product G{zi .,..., Zk+i)A{zi .,..., Z 2 k) is zero. For this purpose, since 

G{zi .,..., Zk+i) is analytic around zero, it is enough to show that there is no 

monomial term in the expansion of A{zi ,..., Z 2 k) of the form 0 ^= 1 ^ 

with cj(i) — 1 < k — 1 for i = 1,..., /c + 1. Since a is a permutation this is 

clearly the case. □ 


4.4.2. A lemma on geometric sums. Let V = diag(i;i,... ,Ud) be a diagonal 
dxd matrix, and M a large parameter. We want to compute the asymptotic 
behaviour of 

d 


(4.37) 


trSym^exp(^F) = 


ki + ...+ka=M 
ki,...,ka>0 


i=l 


This is the coefficient of in the power series expansion of 


1 


det(/-xexp(^F)) ^ = 


1 


^ ' " n?=i(l - 


Lemma 4.12. As M 


oo, 


1 


trSym^^^ exp(— V) = M' 


d-l 


e^^mdxi...dxd + 0{M‘^-^) 


Xl + ...+Xa = l 
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Proof. Dividing by ^ we get a Riemann sum 


Md-i 

k]_+...+kd=M 

ki,...,kd>0 


= JJ ...dxd + 0{^) 

Xl + ...+Xd = l 

Xj >0 


□ 


4.4.3. Example: leading coefficient in the range 0 < c < 1. The leading 
coefficient in Ik{cN,N) (i.e. the coefficient of when 0 < c < 1, can 

be obtained from (14.31) . Thus for 0 < c < 1, we have 7 a;(c) = . We 

now verify that the complicated expression that we got in this section for 
the leading coefficient 7 fc(c), agrees with the above. Note that because of 
the functional equation, Lemma l4.ll we can conclude that this holds also in 
the range (fc — 1)A^ < c < kN 

The leading coefficient of Ik{cN, N) when 0 < c < 1 is 7 a;(c) = 
where 

( 4 . 38 ) s,Ac) = / ■ ■ ■ / 

2k , 

l<2<J<fc j = l 3 

or 

k-\-l<i<j<2k 

2k , 

A{vi,VkfA{vk+i ,..., V 2 k)‘^ n “T 

j=i 

By the residue theorem, in order to compute gkfl{c) we need to find the 
coefficient of 0^=1 ™ ^he expansion of 

+ , VkfA{vk+l, . . . , V2kf- 

Consider the Vandermonde determinant expansion: 

k 

(4.39) A{vi,... ,Vkf = Sgn(cj)Sgn(cj') 

a,a'GSk 

We are looking for terms of the form with a{i) +a'{i) — 

2 < k — 1 for all 1 < z < fc. Since a{i) and cr'{i) are permutations, 
the only such term is the same way, the only possible 

contribution from A{vk+i, ■ ■ ■ ,V 2 k)‘^ to the integral comes from the term 
That means that the term can con¬ 

tribute only a constant. Therefore, the calculation comes down to verifying 
that Jo{cv) = when T = 0. This is indeed the case, since when i; = 0, 
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Jo{cv) is the volume of a — 1 dimensional simplex, which is as 

required. 


5. Justification of Conjecture o 

Our final goal is to sketch briefly a jnstification for Coniecture 11.11 without 
reference to the function field results in the body of the paper. In addition, 
we indicate how to generate a conjecture for the lower order terms in the 
asymptotic expansion as noted at the end of Section [L31 

We start by defining 


(5.1) Qk{a,T) 


la 


r(iogr)^" 


^^2 ^ logT 


+ ii)"C(i + 


la 

logT 




We have the Riemann-Stieljes integral identity. 


OO 

C^(l/2 + ia/logr + R) = j 

— OO 


g-iax/ log (^^xy 


Substituting this into (j5.ip and swapping the order of integration, we find 
that 

OO 

— OO 

Hence, by Fonrier inversion, on average 


(5.2) 


n~<l—v 

(logr)^^-i ’ 


-Ill—1 \ 


Qk{7r/3,T)e^^^^^d/3 


Conjecture II.II now follows from a conjecture of Kosters [26] : 

(5.3) lim Qk{a,T) = ak lim Wk{a,N), 

T—>-oo N^oo 

where we write 

Wk{a, N) = j^ J det(l - e-*"/^t/)^ det(l - dU 

U(N) 

and Ofc is given by (II.9p . (This is a matter of coupling equation (1.2) and 
Conjecture 1.2 of [26] •) We then have, using (14.2p to expand the random 











32 


J.P. KEATING, B.. RODGERS, E. RODITTY-GERSHON AND Z. RUDNIGK 


matrix integral, 

Wk{a,N) = j^ 4(m;Ar)e-2--/^ 

0<m<kN 

0<m<kN 

m>0 

k 

(5.4) ~ y e“^*"“7fc(n) dtt, 

0 


as the last sum is a Riemann sum. 

Setting X = r“ and H = r““i in ()5.2p implies that for x ~ X, on average 


Ai(^; H) 


~ aklk{u). 


We may impose H = by setting u = 1/(1 — 6). The restriction that 
u E [0, k] becomes 5 G [0,1 — 1/A:] and the Conjecture follows. 

The expression (|5.3p follows from conjectures in [7j which relate Qkio(,T) 
to a combinatorial sum, like that in Theorem 14.31 and to a multiple contour 
integral, like that in Theorem 14.101 which include arithmetic factors [26]. 
Specifically, it follows from a leading-order asymptotic evaluation of the 
multiple contour integral that is similar to the calculation given here in 
Section oi A calculation of lower order terms, as in Section 14.21 of the 
present paper, leads to a polynomial of order — 1 in the variable logX 
for the second moment of A^. 
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